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Abstract 

We propose some extensions of the quark potential model to hybrids, fit them to the lattice data 
and use them for the purpose of calculating the masses, root mean square radii and wave functions 
at the origin of the conventional and hybrid charmonium mesons. We treat the ground and excited 
giuonic field between a quark and an antiquark as in the Born-Oppenheimer expansion, and use the 
shooting method to numerically solve the required Schrodinger equation for the radial wave functions; 
from these wave functions we calculate the mesonic properties. For masses we also check through 
a Crank Nichelson discretization. For hybrid charmonium mesons, we consider the exotic quantum 
number states with J PC = 0" 1 , 1 _+ and 2" 1 . We also compare our results with the experimentally 
observed masses and theoretically predicted results of the other models. Our results have implications 
for scalar form factors, energy shifts, magnetic polarizabilities, decay constants, decay widths and 
differential cross sections of conventional and hybrid mesons. 

I. Introduction 

A number of hadron properties are well described by the quark model where mesons have quantum 
numbers J = L © S , P — (— and C = (— 1) L+S , L and S being the quantum numbers for the 

quark- antiquark orbital angular momentum and their net spin angular momentum respectively. The 
states with J PC = H , 1 2 H (for the lowest lying hybrids in the flux tube model) can not be formed 
from a qq pair and hence are not allowed in the quark model. These states are signals for exotic mesons 
(hybrids, glueballs, etc). Quantum Chromodynamics (QCD), describing the interaction between quarks 
and giuonic field, predicts the existence of hybrid mesons containing the excited giuonic field. Thus for 
understanding of QCD, we need to find experimentally testable predictions of the theory for situations 
in which the giuonic field between a quark and antiquark is in an excited state. Thus hybrids are an 
important source of information related to confining properties of QCD, and checking for existence of 
hybrid mesons is very important objective of particle physics. Reviews of the spectrum of excited giuonic 
states can be found in ref. pQ. Recently, a resonance is observed at COMPASS [2] with J PC = 1"+. 
Some other groups like VES [3], E852 [4], and the Crystal Barrel collaboration [5] also observed these 
states. 

Using the Born-Oppenheimer approach, the work of finding implications of QCD for a meson can 
be split into first using the numerical lattice simulations of QCD to fit parameters of a quark antiquark 
potential and then using this potential to calculate its dynamical implications. Even a numerically defined 
potential can be used in this scheme, as in refs. [6] [?]■ These works use some potentials to calculate 
hybrid masses and few radial probability graphs (but presents no qq wave function expressions or uses.). 
What we add to this work is that now we suggest a number of analytical expressions for the excited state 
giuonic field potential between a quark and antiquark and fit their parameters to the lattice data for 
the ground and excited state giuonic field energy values available in ref. [6] for discrete quark antiquark 
separations. For each case, we report a dimcnsionlcss chisquarc and other measures directly telling how 
much fractional error each model has in fitting. We use this variety of potentials to indicate suitable 
ones for different applications, for example for analytical calculations of expectation values in few-body 
wave functions where a numerical approach may have convergence problems. Our full list of chisquare 
and other error measures tells how much extra error is generated in preferring the flux tube model and 
other integrable forms over the ones that can be only numerically used but fit better. After using these 
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potentials to find the quark-antiquark wave functions, we also calculate a number of conventional and 
hybrid meson rms radii and square of radial wave functions at the origin. We continue till writing some 
implications of these for scalar form factors, energy shifts, magnetic polarizabilities, decay constants, 
decay widths and differential cross sections of conventional and hybrid mesons. All this is addition to our 
own results for meson masses and radial probabilities that are in agreement to those reported in ref. [5J 
for the height of the single peak of the probability graph. 

In the present paper, we apply our techniques to charmonium mesons. An advantage of considering 
them is indicated in ref. [8] as 

" The best systems for a hybrid search may be cc or bb where there is a large gap between 
the lowest states and the DD and BB threshod respectively." 

To find the wave functions of conventional charmonium mesons, we use the realistic columbic plus 
linear potential model to solve the Schrodinger equation numerically by using the corresponding quantum 
numbers of mesons. To study hybrids, we repeat the numerical work with the models of the gluonic 
excitation energy mentioned in section III. From the numerically found wave functions, we calculate 
the root mean square radii. These radii can be used to find scalar form factors [5] for charmonium 
mesons, along with energy shifts [10] and magnetic polarizabilities [10] . Thus we have reported some 
predictions about these quantities for conventional and hybrid charmonium mesons. We have also found 
the numerical values of square of radial wave functions at the origin (|i?(0)| 2 ), which can be used to 
calculate the decay constants QT], decay rates [11], and differential cross sections [12] for quarkonium 
states. The predictions about these quantities are also reported for conventional and hybrid charmonium 
mesons. 

In the section II below, we write the Hamiltonian for the conventional mesons. Then we describe the 
shooting method-based numerical procedure to find the solution of the radial Schrodinger equation for 
conventional charmonium mesons. The expressions to find masses, root mean square radii, and squares 
of radial wave functions at the origin (|i?(0)| 2 ) of conventional charmonium mesons are also written in 
this section. In section III, the Hamiltonian is written for hybrid mesons, and then we redo all the 
numerical work as done in section II for hybrids now. The \ 2 an d other error measures for different forms 
of the potential difference between ground and excited state are also written in section III. Results for 
the masses, root mean square radii and |i?(0)| 2 of conventional and hybrid mesons are reported in section 
IV for systems composed of charm quarks and antiquarks. Based on these results, we also include some 
results related to experimentally measurable quantities. 



II. Characteristics of Conventional charmonium mesons 

The Potential Model for Conventional Charmonium Mesons 

In the potential models, the confining potential for QQ system in the the ground state gluonic field 
is mostly used in the form of 

— 4a s , 

-^ + 6r + c, (1) 

with inter-quark distance r. Here, —4/3 is due to the colour factor, a s is the quark-gluon coupling, b is 
the string tension and c is a constant. In above equation, the first term is due to one gluon exchange 
and the second term is the linear confining potential [13] . This potential form provides a good fit to the 
lattice simulations of refs. [131 US]. By including the Gaussian-smeared hyperfine interaction [TB] and 
orbital angular momentum (or centrifugal) term, the potential of the QQ system for the ground state 
gluonic field have following form 

where S c .Sc — ' S ^' S 9 +1 ^ — |, /i is the reduced mass of the quark and antiquark, m c is the mass of the charm 
quark, and S is the total spin quantum number of the meson. For the cc mesons, the parameters a s , b, 
a, and m c are taken to be 0.5461, 0.1425GeV 2 , 1.0946 GeV and 1.4796 GeV respectively as in ref. [TB]. 
The quantum numbers for the conventional charmonium mesons we choose for our study are reported 
below in Table 5. 
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Figure 1: The radial parts of rj,J/ip,h c , and \c meson wave functions as functions of r. Thin solid line 
represents the r] wave function, dashed line represents the J/tp wave function, thick solid line represents 
the h c wave function and points represent the \c wave function. 



Wave Functions and Radii of Conventional Charmonium Mesons 

A conventional meson can be described by the wave function of the bound quark- antiquark state 
which satisfies the Schrodinger equation with potential of eq.(|2]). Radial Schrodinger equation with wave 
function U(r) = rR(r) is written (in natural units) as 

V 2 U(r) + 2^(E -V(r))U{r) = 0. (3) 

Here R(r) is the radial wave function, r is the interquark distance, E is the sum of kinetic and potential 
of quark- anti-quark system, and V(r) and \i are defined above through eq.([2|). 

In quark-antiquark bound state, the wave function must satisfy the boundary conditions U(Q) = 
and U(oo) = 0. For the numerical solution of the Schrodinger equation with the potential of eq.Q), we 
repeatedly generated energy E from -2 to 2 GcV in steps of 0.1 GeV. For each such trial initial energy, we 
used the Newton method [17] to select, if any, the energy for which the numerical solution of Schrodinger 
equation became zero at infinity. To obtain these numerical solutions, we used the RK method [18j 
with using any arbitrary integer value of U'(0). For different values of U'(0), normalized solutions of the 

Schrodinger equation, obtained by multiplying the solution with the normalization constant ( , ), 

y/JU 2 (r)dr 

remain the same. These energy eigenvalues plus constituent quark masses are taken to be the cc mesons 
masses (in natural units). It is found that our results for conventional charmonium mesons agree with 
the Table 1 of ref. 16.. This supports the reliability of our method. We also checked the consistency 
of our method by 1) getting a 100% overlap of our HU and EU and 2) by calculating the masses of 
conventional mesons by the Crank Nichelson Discritization and finding that masses obtained by both 
of the methods are identical. The Fig.l shows the dependence of 77, J/if), h c , and \ c normalized radial 
wave functions on the radial separation r. The quantum numbers (L and S) for these particles are given 
below in Table 5. These graphs show that the radial wave functions of 77, h c , Xa & n d J/tp have the same 
properties as that of hydrogen atom radial wave functions, i.e. they behave as r L for small inter quark 
distances and decrease exponentially at large inter quark distances. Thin solid and Dashed lines graphs 
are for L = 0, so these graphs are similar to r°exp(— r). Thick solid and points are for L = 1, so these 
graphs are similar to r exp(-r). As L increases, the wave function's peak goes away from the origin. 
This means that centrifugal term has more effects on wave function than that of the hyperfine term. One 
possible reason is that we are dealing with heavy quarks so the l/wi^ factor (shown in eq.©) of the 
hyperfine term becomes very small. 

The normalized wave functions are used in the further calculations for root mean square radii and 
square of radial wave functions at origin. To find the root mean square radii of the cc mesons, we used 
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Figure 2: Graphs of potential energy differences between ground and excited state. Points represent 
the data taken from ref. [5J, thick solid line represents potential difference model — , thin solid line 
represents the potential energy difference ^ + Ax exp~ Br , dashed line represents the potential difference 
A x exp~ Br ° 1897 , and squared points represents the potential energy difference A x exp~ Br ' 2 . 



the following relation: 



(r 2 ) = U*r 2 Udr. (4) 

In ref. |19j . for normalized wave function 

U'(0) = R(0) = n/4^V(0) (5) 

is used and we use this prescription. Thus the derivative of U(r) at r — is calculated to find |i?(0)| 2 . 
|i?(0)| 2 is used in many applications of high energy physics as mentioned in section I. 



III. Characteristics of Hybrid Charmonium mesons 

The Potential Model for Hybrid Charmonium Mesons 

The centrifugal factor for the hybrid mesons is written in refs. [2"Ul [2Tj as 

L(L + l)-2A 2 +< J 2 > 

2^ (6) 

where A is the projection of the total angular momentum J g of the gluonic field. The states with 
A = 0, f , 2, 3, , .... are usually represented by the capital greek letters E, II, A, ... respectively. We are 
interested in finding the masses and root mean square radii of the hybrid states H , 1 ^ and 2 H . 
These states can be generated from the II U potential. For the fl u potential, < J 2 >= 2 and A = 1 [21] . 
Therefore — 2A 2 + (J 2 ) = 0, so centrifugal factor for the hybrid mesons is L(L + l)/2mr 2 . In ref. [21] 
J = L © S, P = e(-l) L+A+1 , and C = eri(-l) L+A+s with e,r) = ±1. Therefore with same quantum 
numbers (L,S), different J PC states are possible. L and S for these hybrid J states are shown in Table 
6 (as given in ref. 22 ). For the excited part of quark antiquark potential ir/r is used in the flux tube 
model [23 . This form of excitation energy is only valid at large inter-quark distances. In comparison, we 
suggest and evaluate excited potential energy in forms which are valid for smaller distances as well. For 
this purpose, we get the potential energy differences (e,) between ground and excited states for different 
r, values from the lattice simulation reported in Fig. 3 of ref. [BJ, and calculate % 2 with a variety (ir/r, 
Axexp ,Axexp ,^+Axexp , and - -I- A x exp ) or ansatz by fitting parameters 
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Table 1: Our calculated x 2 for the data of potential difference to our suggested models {V g (r)) with best 
fit parameter's values. 













Parameters 




ansatz 


x 2 


X 2 /D 


X 


A 


B 


c 


7 






GeV 2 




GeV 


GeV 






ir/r 


0.2305 


0.3268 


0.5285 










A x cxp(-Br 7 ) 


0.0857 


0.1215 


0.2773 


1.8139 


0.0657 




2 


A x cxp(-Br 7 ) 


0.0004 


0.0005 


0.0205 


17.9325 


2.5195 




0.1897 


f + A x cxp{-Br~<) 


0.00f2 


0.0017 


0.0331 


1.2448 


0.1771 


0.3583 


1 


l+A x exp(-Br^) 


0.0003 


0.0004 


0.0132 


3.4693 


1.0110 


0.1745 


0.3723 



appearing in each ansatz. Dimensionless x 2 is defined as 



with i = 1, 2, 3, n being number of data points. Here V g (r) is a model of the potential energy difference 
between the ground and excited state. We tried 

V g( r ) = 7r / r > 

V g (r) = A x exp(-Br 2 ), 

V g (r) = Ax cxp(-Br^), 

V g (r) — - + A x exp(—Br), and 

V g (r) = I + A x exp(-BrT). 

The parameters A, B, c, and 7 are found by fitting these models (V g (r)) with the potential difference 
data of ref. [6], and reported in Table 1. \ 2 , x 2 /D an d X are also reported in Tables 1 with X 2 /D = 



X 2 /(Number of data points) and x — 

It is noted that the x 2 for ir/r is greater than all other potential difference forms. Fig. 2, representing 
the graphs for different forms of potential difference, also shows the same behaviour. The x 2 of the 
gaussian gluonic potential (A x exp~ Br ) is less than but larger than all other forms. A x exp~ Br 
is a smeared form of con * > as written in appendix of ref. [24]. The potential difference in this form 
(A x exp~ Br2 ) has an advantage that it can be easily used in dynamical applications. For example, 
the expectation value of this part of the potential energy in a Gaussian wave function of the quadratic 
confining potential is given in terms of error function differences even if one uses antiderivatives to evaluate 
the definite integrals in it. Or, for usual infinite limits it simply multiplies in the integrand with already 
Gaussian meson wave functions to keep the integrand as Gaussian, whose well known integral can be 
written by inspection. (The expression for the expectation value can be minimized with respect to the 
chosen parameters of the wave function to find the ground state energy and wave function using the 
variational method.) This calculational advantage may be trivial for a two-body problem. But if one 
has to evaluate an expectation value for few or many body problem (or for a minor variant used in the 
resonating group method based treatment [25] [26] of a system of two quark and two antiquarks), we 
have to evaluate an integral of a high order whose direct numerical evaluation may have convergence 
problems as in ref. |26j and a Gaussian integration by inspection may well be the only practical option. 
The need to keep the multi-dimensional integrals as Gaussian, giving importance to the A x exp~ Br 
form, becomes even more prominent when when wave functions of the conventional mesons are replaced 
for the respective problems by those of the hybrids as in ref. [25] . 

The analytic Gaussian expectation value of A x exp(— Br) term is similarly given in terms of error 
function differences. (Or, its product in integrand with Gaussian wave functions can be converted to a 
new Gaussian integrand using a completing of square.) The x 2 fo r the potential difference in form of 
- + A x exp(— Br 3723 ) is very less. The - term in this form can be used for analytical expressions 
of expectation values in the above mentioned Gaussian quark antiquark wave functions of quadratic 
potential, resulting in differences of the exponential integral functions for the most analytical way of 
finding the expectation values and integrals of resonating group method [25] |26j . But a similar fully 
analytical route for Gaussian expectation values and resonating group integrals of the A x exp(— _Br 7 ), 
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Table 2: Our calculated \ 2 for the first excited state data of ref. [6] to the model used in ref. [7J with 
best fit parameter's values. 



Excited Potential 


x 2 


X 2 /D 


X 


bo 


Par am 
bi 


eters 
b 2 


co 






GeV 2 




GeV z 


GcV a 


GcV 4 


GcV 


co+ 

Vbo + b\r + b 2 r 2 


0.00096 


0.0030 


0.0226 


58.00f6 


4.4896 x 
10 -io 


0.2859 


-6.1814 




Figure 3: Graphs of excited state potential energy in the form of Co + V&o + b\r + b 2 r 2 along with data 
of ref. [6]. Solid line is for the potential cq + V^o + b\r + b 2 r 2 and dots are for data of ref. [6]. 



for 7 = non-integer number, is not available, and this can lead to convergence problems [26] when we 
integrate numerically integrals of high dimensions. 

In ref. [7J, [21] excited state potential (not the difference) is used in the form of 



Cfl 



+ \/bo + hr + b 2 r 2 . (8) 



By fitting this potential with data of first excited potential taken from Fig. 3 of ref. [6], we calculated 
the parameters bo, b\, b 2 , c . We also calculated x 2 , \ and x 2 /D an d reported in Table 2 along with 
the parameter values. The fit of the data with this excited state potential is shown in our Fig. 3. As for 
the the analytical calculations (for finding expectation values etc.), these are also not possible with this 
potential form and thus for many applications it has to be replaced by others of higher chisquare. 
Fig. 4 shows the behaviour of string potential [37J 



a 2 r 2 + 2na{N +'^) + w q (9) 

with N = 0,1, 2. In this figure, the points represent the data taken from ref. (ij . The parameter's values, 
X 2 and other error measures of the string potential with the data (of excited potential) of ref. [6] is 
reported in Table 3. The parameter's values are calculated by fitting the data of excited potential with 
the string potential. The analytical calculations (for finding expectation values) are not possible with 
this string potential form as well. 

Wave Functions and Radii of Hybrid Charmonium Mesons 

Now, we can write the quark antiquark potential in excited state gluonic field as 

m _ zigi + br + "g- ( » ^ w ^ + " < J i > + w . do, 

3r 9mj ^/tt 2[ir z 
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Figure 4: Graphs of potential energy in the form of a 2 r 2 + 2ttcj{N + |) + w q . The solid line is for 

N = and small points are for N — data of ref . [6] . Thick solid line is for N = 1 and medium size 
points arc for N = 1 data of ref. [B] . Dashed line is for N = 2 and large size points are for N = 2 data 
of ref. |B]. 



Table 3: Our calculated \ 2 with best fit parameter's values. 



Excited Potential 


x 2 


X 2 /D 


X 


N 


Parameters 
w q 


a 






GeV z 


GeV 




GeV^ 


GeV 4 




0.5536 


0.6633 


0.6932 





5.9884 x 10" 9 


0.0601 


^Ja 2 r 2 + 2ttct(N + §) + w q 


0.0112 


0.0348 


0.0899 


1 


3.5044 x 10~ 7 


0.1629 




0.0190 


0.0740 


0.1259 


2 


1.0531 x 10" 5 


0.1627 
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Table 4: Our best fit parameter's values and overlaps of numerically solved wave functions to the wave 
function form written in eq. (jll|) . 







L = 1 and 5=1 






L = 2 and 5=1 






Parameters of excited state 


Overlap 


Parameters of excited state 


Overlap 




function written in eq.(|ll[) 




function written in eq.pTj) 




ans'dtz 


n 


P 




n 


P 




7r/r 


0.0047 


0.0556 


0.9997 


0.0032 


0.0451 


0.9989 


A x cxp(— Br 2 ) 


0.0031 


0.0439 


0.9964 


0.0025 


0.0391 


0.9905 


A x cxp(-5r - 1897 ) 


0.0093 


0.0822 


0.9903 


0.0048 


0.0565 


0.9998 


f + A x exp(-Br) 


0.0070 


0.0698 


0.9899 


0.0036 


0.0482 


0.9999 


f + A x cxp(- J Br ' 3723 ) 


0.0087 


0.0794 


0.9903 


0.0045 


0.0547 


0.9998 


Co + V^o + b\r + b 2 r 2 


0.0163 


0.1135 


0.9921 


0.0090 


0.0807 


0.9999 



V g (r) is denned above after eq.Q. 

Using this excited state potential of eq. (fT0f along with the above mentioned values (after eq.©) of A 
and < J 2 >, the energy eigenvalues and the corresponding wave functions are found by using the same 
technique as employed for conventional mesons (mentioned in section II). As before, these eigenvalues plus 
constituent quark antiquark masses are taken to be the masses of hybrid mesons. Then we normalized the 
wave functions and found the root mean square radii of hybrid mesons by using cq.Q. The normalized 
radial wave functions for charmonium hybrid mesons are graphically represented in Fig. 5 and Fig. 6. 
The overlaps of our numerically calculated radial wave functions (U = rR) for the excited states and a 
modified gaussian wave function ansatz 

ip = n r 2 exp(— pr 2 ) (11) 

multiplied by \/4irr are written in Table 4 in such a way that U = rR = y/farrt/). 
The normalization of the gaussian wave function gives 

n = (423p3)(1557T3). (12) 

The numerical value of p is found by fitting this function with the data of numerically calculated wave 
function, and written in Table 4. 

The Fig. 5 and Fig. 6 show the wave function dependence on L and 5. Therefore the masses and root 
mean square radii of H ,1 h and 2 H J PC states also depend on the quantum numbers L and 5. 
|i?(0)| 2 is found for hybrid mesons using eq.flS). The Fig. 5 and Fig. 6 also show that the peaks of the 
graphs for radial probability density are in agreement with the peak of radial probability curve drawn in 
ref. 0. 

IV. Results and Conclusions 

1- For conventional mesons, our calculated masses and root mean square radii are reported in Table 
5 along with the experimental and theoretical predictions of the other works. We observed that our 
results are in good agreement with the experimental and existing theoretically predicted values, which 
shows the validity of our method. Quantum mechanically, when L increases, centrifugal barrier increases 
so particles become less bound implying increased root mean square radii. Our calculated root mean 
square radii are in agreement with this expectation. 

2- With the parameters (given in Table 6) for the 0+",l -+ and 2+ - J PC states, masses and root 
mean square radii are calculated for the charmonium hybrid mesons. In Table 6, masses are calculated 
using the excited state gluonic field potential in the above mentioned forms. For comparison with earlier 
works, masses of cc hybrid mesons with H , 1 h and 2 H J PC states are given in Table 7. In Table 
8, root mean square radii are calculated by taking the excited state potential in the coulomb plus linear 
plus additional excited potential. In Table 9, masses and root mean square radii are reported for the 
excited potential in the form of cq + -\/&o + b\r + b 2 r 2 . 
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Figure 5: Hybrid charmonium meson radial wave functions for H , 1 h and 2 H J PC states for L=l 
and S=l. The wave function with potential in the form of coulombic plus linear plus Aexp{— Br 2 ) is 
represented by the solid line. Wave function with coulombic plus linear plus ir/r potential is represented by 
square points. Wave function with coulombic plus linear plus A x Exp(—B r 1896 ) potential is represented 
by dashed line. Wave function with coulombic plus linear plus ^ + Ax Exp(—B r) potential is represented 
by thin solid line. Wave function with coulombic plus linear plus f + A x Exp(—Br 03723 ^ potential is 
represented by thick dashed line, and the wave function with excited potential in the form of cq + 
Vbo + b\r + b^r 2 is represented by points with lines. 



Table 5: The experimental and theoretical masses and theoretical root mean square radii of some con- 
ventional charmonium mesons. The experimental mass is the average PDG |16j and rounded to 0.001 
GeV. Our calculated masses are rounded to 0.0001 GeV. 



Meson 


L 


S 


Our calculated 


Theor. mass [16] 


Exp. mass 


our calculated 


Theor. ^(r 2 ) [28] 








mass 


with NR potential 




VF) 


with 










model 






potential model 








GeV 


GeV 


GeV 


fermi 


fermi 


Vc 








2.9816 


2.982 


2.9792 ±0.0013 [16 


0.365 


0.388 


J/1> 





1 


3.0900 


3.090 


3. 09687 ±0.00004 [16] 


0.414 


0.404 


h c 


1 





3.5156 


3.516 


3.525 ± 0.00055 [16 


0.674 


0.602 


Xc 


1 


1 


3.5246 


3.556 


3.55618 ± 0.00013 29^ 


0.685 


0.606 
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Figure 6: The same figure as Fig. 4 but with L = 2 and S = 1. 



Table 6: Our calculated masses of cc hybrid meson 0+-,l- + and 2+~ J PC states. 



jPC 


L 




A 


<Ji> 


Excited potential as coulombic plus linear plus 












n/r 


Ax 
exp(— Br 2 ) 


Ax 

exp(-Br - 1897 ) 


^ + Ax 

r 

exp(-Br) 


^ + Ax 
exp(-£r°- 3723 ) 












GeV 


GeV 


GeV 


GeV 


GeV 


0+-,l-+,2+- 


1 


1 


1 


2 


4.3571 


4.0619 


4.2680 


4.2733 


4.2694 


1-+.2+- 


2 


1 


1 


2 


4.4632 


4.1433 


4.4632 


4.4258 


4.40796 



Table 7: The mass predictions of 1 + , + and 2 + states of other works. 



Predicted masses (GeV) 


models 


1-+ 


o+- 


2+- 




w 3.9 [30 a 






bag model 


4.2-4.5 31 






flux tube model 


4.19 ± sys. error 32 [33 
4.7 [3U 


«4.5 H 
4.58 [32] 




heavy quark LGT 


4.1-4.5 






QCD sum rules 


4.369-4.420 [35, 36, 37J 


4.714(52) [37 


4.895(88) 38, 


quenched lattice QCD 



10 



Table 8: Our calculated root mean square radii of cc hybrid meson + , 1 + and 2+ J PC states. 



J PC 


L 




A 


< J 2 > 


(r 2 ) with excited potential as coulombic plus linear plus 












n/r 


Ax 
exp(— Br 2 ) 


Ax 

exp(-Br - 1897 ) 


£ + Ax 

r 

exp(-Br) 


^ + Ax 
exp(-Br°- 3723 ) 












fm 


fm 


fm 


fm 


fm 


0+-,l-+ 2+" 
l-+,2+~ 


1 

2 


1 
1 


1 
1 


2 
2 


1.1061 
1.2280 


1.2458 
1.3203 


0.9110 
1.0988 


0.9881 
1.1883 


0.9272 
1.1160 



Table 9: Our calculated masses and root mean square radii with excited potential Co + V&o + b\r + b 2 r 2 
of cc hybrid meson with H , 1 h , 2 H J PC states. 
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<J 2 9 > 


masses 


V(r 2 ) 












GeV 


fm 


0+-,l" + ,2+- 


1 


1 


1 


2 


4.9503 


0.7748 


1-+.2+- 


2 


1 


1 


2 


5.1693 


0.9185 



3-For conventional mesons |i?(0)| 2 is reported in Table 10. Each of |i?(0)| 2 of cc hybrid mesons for 
H , 1~ + and 2 H J PC states is equal to zero by our calculation and this result agrees with ref. 39 which 
writes "models of hybrids typically expect the wave function at the origin to vanish" . We also noted that 
the masses and root mean square radii of the hybrid mesons are greater than ordinary mesons of same 
flavour and quantum numbers (L and S). Since H , 1 ^ 2 H states are not possible with quark model 
quantum numbers, so we can not compare these J PC states with conventional mesons. 

In ref. [9], the scalar form factor is written as 

r.(<) = i+ l -(r 2 y s t + o(t 2 ). (13) 

6 

In ref. |10j . energy shift and magnetic polarizability are written as 

AE n = {n \ —L 3 + — + ■ — \n) + S„ . (14) 

8(i \4:fx (Mi + [i 2 ) 32 E n - E n , 

Here the symbol L 3l H, e, m are used for the angular momentum, magnetic field, charge and mass of the 
quark, Jl = ^, 8 is the angle between H and relative co-ordinate r, q = ei — e%, and e = e\ + ei. In above 

equation, the term having (n \ H r 3 " 9 \ n!) is related to square of root mean square radii. 

P = -^(i + ^—]{r 2 ). (15) 
24 V4^ (ii+mj 

In above eas. (|13til5|) the root mean square radii is in the numerator, therefore we predict that magnitudes 
of scalar form factor [9], energy shift [ID] , and magnetic polarizability 10 for hybrids are greater than 
those for conventional mesons of the same quantum numbers {L and S) . 

By parametrizing the excited state wave function written above in eq. (|ll[) . we get 



(r 2 )=Jrr^dr= - 28 y f i- . (16) 

Here n is a function of p, as given by eq. (fl"2"]) . Substituting the result of eq. flTBl in eqs. (|13U5[) . the scalar 
form factor and magnetic polarizability become 

- , , 1 15n 2 VW2 , o N 

r„(t) = i + - 10Q V 7/ ; t + (t 2 ), (17) 



6 128p 7 /2 
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Table 10: |i?(0)| 2 of cc meson 



Meson 


L 


S 


our calculated normalized \R(0) 2 
(fm)- 3 


V 








2.5910 







1 


2.1789 


h c 


1 





< io- lu 


Xc 


1 


1 


< io- n 



1 (e 2 g 2 \ 15n 2 VV2 
P 24^4/7 128pV2 ■ ^ 

Numerically calculated values of n and p are written in Table 4 for different forms of excited state 
potentials for different L and 5. 

As we mentioned above, |i?(0)| 2 is equal to zero for hybrid mesons. Using this result, we can predict 
that decay constants [11] , decay rates [11] , and differential cross sections [12] of hybrid mesons are zero 
as these quantities are proportional to |i?(0)| 2 as written in these references. 
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